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Abstract 

We develop a new integration technique allowing one to construct a rich manifold of 
particular solutions to multidimensional generalizations of classical C- and S*-integrable 
Partial Differential Equations (PDEs). Generalizations of (l+l)-dimensional C-integrable 
' and (2+l)-dimensional .S-integrable iV-wave equations are derived among examples. Ex- 
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amples of multidimensional second order PDEs are represented as well. 



1 Introduction 

After the original work pQ, the integrability technique was intensively developing. At present 
time, it covers large class of nonlinear Partial Differential Equations (PDEs) applicable in 
different branches of physics and mathematics. One should mention a few most popular in- 
^-j- ! tegration methods, such as (a) the linearization by some direct substitution, for instance, 
by the Hopf-Cole substitution [2J and by its multidimensional generalization [3] (appropriate 
nonlinear PDEs are referred to as C-integrable PDEs, [U El El El El E]); (b) the method of 
OS characteristics JTU] and its matrix generalization [HI [12] {Ch- integrable PDEs); (c) the inverse 
spectral transform method [U [131 EH] > t ne dressing method [151 DEI El [HI EH] an d Sato ap- 
proach [20] (S'-integrable PDEs). S- integrable nonlinear PDEs are most applicable in physics. 
We recall a few types of these equations: the soliton equations in (l+l)-dimensions, such as 
the Korteweg-de Vries (KdV) [U EI] and the Nonlinear Shrodinger (NLS) [22] equations; the 
soliton (2+l)-dimensional equations, such as the Kadomtsev-Petviashvili (KP) [23J and the 
^ Davey-Stewartson (DS) [21] equations; the self-dual type PDEs having instanton solutions, 
like the Self-dual Yang-Mills equation (SDYM); PDEs assotiated with commuting vector fields 
[23 EHl E3 EH EH EQl El]. Nevertheless, the class of completely integrable nonlinear PDEs 
is very restrictive. Thus, extensions of the integrability technique on new types of nonlinear 
PDEs is an actual problem. 

In this paper we suggest a new version of the dressing method allowing one to construct a rich 
manifold of particular solutions to new class of nonlinear PDEs in any dimension. The novelty 
of this class of PDEs is that integrability technique does not generate commuting flows to them 
in usual sense, unlike all methods mentioned above, where any nonlinear PDE appears together 
with the commuting hierarchy of nonlinear PDEs. We show that our algorithm may provide 
arbitrary functions of m — 1 independent variables in the solution space to m-dimensional PDEs, 
which suggests us to consider these PDEs as candidates for complete integrable PDEs. However 
we do not represent rigorous justification of complete integrability. 

To anticipate, we give simple examples of nonlinear PDEs for the matrix field V, derived in 
this paper. 
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1. The system of first order D-dimensional PDEs, 



D 



J2{Vt m + VC {m) V)B^ =0, (1) 



m=l 



where B^ and are some constant matrices. This is multidimensional generalization of 
C-integrable (l+l)-dimensional nonlinear iV-wave equation. 

Simple example of this equation corresponds to D = 2, t± — x, t 2 = y, B% = B^ = 0, 



Then eq.(JTJ reads 



u x + u 2 C^ ] +pqC { 2 ] = 0, (3) 
p x +p(uC^ ] +vC { 2 l) ) = 0, 
v y + v 2 C {2) +pqC (2) = 0, 
q y + q(uC[ 2) + vCi 2) ) = 0, 

which reduces to the Liouville equation 

/ a , = cf cfexp(2/), (4) 

if = = 0, q = p = ef. Here and below, B^ and Ca mean the ath diagonal elements 
of the diagonal matrices I?( m ) and C ( - m * ) respectively. 
2. The first order .DiZ^-dimensional PDEs 

L(mi) ( V ^m 2 + VC^^V) = 0, (5) 

rri2=l mi=l 

where C^ m ^ and are some constant matrices. The simple example corresponds 

to Di = D 2 = 2, Li" 11 - 1 = R^ 12 ^ = if a ^ m x and (3 ^ m 2 respectively, tu = x, t\ 2 = y, 
t 2 \ = z, t 2 2 = t. Let V be given by eq.(j2J), then eq.(jH]) yields: 

u x + M 2 C; n) + pqC { 2 11] = 0, (6) 
p z +p(uC {2l) +vC {21) ) = 0, 
Vt + v *CW+pqCi 22) = Q, 
q y + q(uC[ 12) + vC ( 2 12) ) = 0, 

which reduces to the following four-dimensional generalization of Liouville equation (j3j) 

9xy = cf>cfV+', u = cf^cf (7) 

if = & 22) = Cf ] = Cf> =0,p = ex P f,q = exp g. 
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In particular, eq.flH} reduces to the following D (D — l)/2-dimensional PDE, D\ = D 2 = D 

D Do 

£ E [i {L {m) V T ^L^ - L^V T ^L^) + (8) 

7712 = 1 m l = ! 

mo > m\ 

^{m^y (j{m 1 m 2 )y j^{m 2 ) _ jjm 2 ) y(j(rn 1 m 2 )y jjpii ) j _ 

where C^ mim2 ^ = — C^" 12 " 11 ^ 1 ; C^ mim2 ^ and L^ mi ^ are constant diagonal matrices and + means 
hermitian conjugate. This equation has a physical mening describing interaction of no(rio — 1)/2 
waves if is no x hq matrix. This is multidimensional generalization of S'-integrable (2+1)- 
dimensional nonlinear iV-wave equation. 

3. The system of second order D-dimensional PDEs, 

D 

( Vt ^m + (VC (n) V) tm + VC^V tn + VC^VC^V) B^ = 0, (9) 

m,n=l 

where B^" 1 ^ are constant diagonal matrices. The scalar version of this equation reads 

D D 

£ V tntm B^ + Vt m C^ + V 3 C = 0, (10) 

m,n=l m=l 

D D D 

n=l m=l m,n=l 

Note that not all constant coefficients may be arbitrary in the above nonlinear PDEs. 
Constructing particular solutions, we will reveal some relations among coefficients. 

Particular examples of nonlinear PDEs, such as eqs. fl3|6|8lfT0|) . will not be considered in 
this paper. Instead of this, we concentrate on the dressing algorithm allowing one to derive 
general equations, such as eqs. (fTTl8ll9l . and study the richness of available solution space for 
them. 

The structure of this paper is following. We will derive generalization of the classical C- 
integrable first and second order nonlinear PDEs in Secj2] with eqs.Q and Q as particular 
examples. Richness of the solution space to eq.flT]) will be discussed briefly and explicite par- 
ticular solutions to this equation with D = 3 will be given. Generalization of the classical 
S'-integrable nonlinear PDEs will be consider in SecJHlwith eqs.(J5]) and (JS} as particular exam- 
ples. Richness of the solution space to eq. <^ will be discussed briefly and explicite particular 
solutions to this equation with D\ = D 2 = 2 will be given. Conclusions will be represented in 
Secgl 

2 Generalization of C-integrable nonlinear PDEs 

2.1 Starting equations 

Our algorithm is based on the following integral equation 

P(jj) * X (j*, A; t) = W(jjl; t) * xGu, A; t) + W(X; t) = W(/x; t) * ( X ( M) A; t) + I x (m, A)) (11) 
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where -P(/i), xCz-tjAjt), W(X;t) are n x n matrix functions of arguments, U are independent 
variables of nonlinear PDEs, t = (ti, . . . , to), rank(P) = no, A, /i, v are complex parameters. 
Here * means the integral operator defined for any two functions /(/i) and g(/i) as follows: 

(12) 

/0) * #0) = J f{v)g{v)d&i{v), (is) 

and fii(/i) is some measure. We also introduce unit Xi(A,/i) and inverse / -1 (A, /x) operators as 
follows: 

/(A, v) * Z^, n) = X x {\, v) * f(v, n) = /(A, /i), (14) 
/(A, r/) * /">, fi) = f-\\, v) * f{v, ft) = Z a (A, /i). 

We introduce parameters ij through the function x(A,/i;t), which is defined as a solution to 
the following system of linear equations: 

XtJX, /x; f) = A^(\, u) * /i; t) + A^(X)P(u) * x (i/, //; 1) + ^(A, /i), m = 1, . . . , £>,(15) 

where A^(A, i/) and v4^ m ^(A) are no x no matrix functions of arguments. An important re- 
quirement to eq.f JTTj) is that it must be uniquely solvable with respect to W(X;t), i.e. operator 
*( X (/x, A; t) + Xi(/i, A)) must be invertable. 

Matrices and A^ m ' may not be arbitrary. They have to provide compatibility of the sys- 
tem (I15p. We will show that there are two different methods which provide this compatibility. 
The first one (Sec J2.2p yields classical C-integrable nonlinear PDEs, linearizable by the mul- 
tidimensional version of the Hopf-Cole transformation [2], while the second method (SecfjO]) 
yields a new type of nonlinear PDEs whose complete integrability is not clarified yet. However, 
our algorithm supplyes, at least, a rich manifold of particular solutions to these PDEs. 

The following theorem is valid for both cases. 

Theorem 2.1. Matrix function W(X;t), obtained as a solution to the integral equation 
ffTTj) with x defined by eq.f lT5|) satisfies the following system of compatible linear equations 

E^(X; t) := W tm (X; t) + V^(t)W(X; t) + (W(n; t) - P(u)) * A™(n, A) = 0, (16) 
m »({) = (%;t)-P(/i))*i (m) (/i), m = l,...,D. (17) 

Proof: To derive eq.tTIB"]). we differentiate eq.f lTTj) with respect to t m . Then, in view of 
eq. (|T5|) . one gets the following equation: 

E {m) (^ t) * ( X (/i, A; t) + Zi( Mj A)) = (18) 

where is defined in eq. flTBl) . Since operator *( X (/i, A; t) +Xi(/i, A)) is invertable, eq. ljTB"]) 
yields: E^ m \n;t) = 0, which coinsides with eq. ljTB"]) . ■ 

Remark: Following the classical integrability theory, we refer to eq. (I16p as the linear equation 
for the function W(X;t). However this is not correct, because functions V {m \t) are defined in 
terms of H^(A; t) by eq. (fT7|) . Thus, strictly speaking, eq. ffTB"]) is a nonlinear equation for W(X] t). 

System ffTB"]) is overdetermined system of compatible linear equations with potentials (t) 
in analogy with the classical integrability theory. In the classical theory, nonlinear PDEs for 
potentials may be obtained as compatibility conditions for the appropriate overdetermined 
linear system. However, this approach does not work in our case because of the last term in 
eqs. ffTB"]) . Instead of this, we suggest a different method of derivation of nonlinear PDEs, see 
SecsJUandO 

Now we analyze two methods that provide the compatibility of system (|T5l) and derive 
nonlinear PDEs assotiated with each of them. 



4 



2.2 First method: classical C-integrable nonlinear PDEs 

In this subsection we write the compatibility condition of eqs. (fT5|) as follows: 

A™ (A, u) + A™ (A)P(i/)) * X t n (u, fi) = (A.™ (A, u) + A™ (A)P(z/)) * X t m (*/, //) , V n, m. (19) 
Substituting eqs. f|T5|) for derivatives of x into eq. f|T9|) we obtain the following equation: 



(L (m) * L (n) - L (n) * L (m) ) * x + ^ M * A {n) - L (n) * A {m) = 0, (20) 
L( m )(A,/i) = A^(A, / u) + i (m) (A)P(Ai), . 

Let eq. (l20|l be satisfied for any function x(\ U>\ t) (which is a solution to the system (TL5]) ). Then 
eq.f l20p is equivalent to two following equations relating matrix functions A^ m \ A^i and P: 

£M * _ lW * = o, (21) 

L M*L<ri-LW*L<ri = eq JP (22) 
(L (m) (A, i/) * A (n \u) - L (n) (A, v) * i (m) (z/))p(^) = 

Since rank(P) = n , eq. (!22j) is equivalent to the following one: 

£M * _ * = o ^ (23) 
* _ * A {m) = A {n) P * A (m) - A (m) p * A (n \ n, m = 1, . . . , D. 

Eqs. (121j) and ( 1231) represent two constraints on the functions A^(A, //) and A^(A). 

Now we have everything for derivation of nonlinear PDEs for the fields IA m )(t). For this 
purpose, let us consider the following combination of eqs. (JIB]): 

E {m \\- 1 ) * A {n \X) - P (n) (A; t) * i {m) (A), (24) 

which yields, in view of eq. (j23~j) 

v£\t) ~ v t\t) + V {m \ t )V {n) (t) - V {n \t)V {m \t) = 0. (25) 

This equation is known to be linearizable by the Hopf-Cole transformation (2j |3]: 

Ef :=^ tn {t) = ^it)V {n \t), (26) 

where ^(t) is an arbitrary n x n matrix function of all variables £j. The presence of an 
arbitrary function ^(t) is assotiated with the fact that the system of nonlinear PDEs (j25p is 
not complete. One needs one more equation relating V^ n \ n — 1, . . . , D. 

To derive this additional equation we introduce either additional relations among 

AW a nd 

in our algorithm based on eq.f fTTj) or an additional linear PDE for \I> in the classical 
algorithm based on eq.f )26|) [3]. 
For instance, let 

D D 

A(m) (A, v) * A {n) {u)B (m) = - A {m) (A)P(i/) * i (n) {u)B (m) (27) 

m=l m=l 
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and/or 



D 



J2% m B^ = 0, (28) 



m=l 



where are no x n arbitrary constant matrices. Then both the combination of eqs. (|T6|) . 

TZ=i( E{n) ) B[m) \ and the appropriate combination of eqs.(|25|), * _1 ELi( £ ?)tm 5(ro) ' y ield 
the same nonlinear equation for V^: 

D 

{yt ] + V {n) V {m) ^j 5 (m) = 0, n = 1, . . . , D. (29) 

m=l 

This A^-wave type equation is supplemented by constraints ( |25l) [3]. 

In particular, introducing reduction V^ m \t) = V(t)C^ (where are no x no constant 
matrices) we reduce the system ( |25i) into the following one: 

V tm C {n) - V tn C {m) + VC {m) VC {n) - VC {n) VC {m) = 0, V n, m, (30) 

which is a (l+l)-dimensional hierarchy of commuting C-integrable iV-wave equations. 

Higher order nonlinear PDEs may be obtained in a similar way introducing appropriate 
equations instead of eq. fl27|) and/or eq.f l28|) . 

C-integrable PDEs will not be considered in this paper. 



2.3 Second method: new class of nonlinear PDEs 

In this section we represent another way to provide the compatibility of eqs.( TT5l) . As a result 
we obtain a new class of nonlinear PDEs together with the rich manifold of particular solutions. 
In particular, there are solutions in the form of rational functions of exponents. 
Let us use the following representation of A^ m \\,fi) and A^ m '(X): 

A( m \X,n) = a( m \X,v)*(3 {m \v,n), A^ m \X) = a {m) {X, v) * P (m) {u), (31) 

where a^ m '(X, v\ p( m '(is, jj) and ^ m >{y) are no x no matrix functions of arguments. Here we 
introduce one more integral operator * defined for any two functions f(fi) and g(fi) as follows: 

f(ji)*g(ji) = J f(ji)g(jj)dn 2 (ji), (32) 

where fi 2 (/x) is some measure. We also introduce unit Z 2 (A,/i) and inverse /~ 1 (A,/i) operators 
as follows: 

/(A, v) * l 2 (u, fi) = X 2 (A, u) * f(u, fi) = /(A, fi), (33) 
/(A, v) * f-\u, /i) = / _1 (A, v) * /(!/, M ) = X 2 (A, /i). 

Now let us write ec[.(H5|) using representation f[3"Tj) in the following form: 

Xt m (X,fi;t) = a( m \X,v)*^ m \v,^,t). (34) 

Here 

e (m) (A, n- 1) = (> m )(A, u) + P {m \X)P(u)) * x (y, /,; t) + /^(A, „), m = 1, . . . , D. (35) 
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It is obvious that the compatibility condition of the system ( Tl5i) is equivalent to the compati- 
bility condition of the system (134|) which reads 

a (m) (A, v) * 4 m) (v, //; t) — (A, v) * $ A*! *) , V n, m (36) 
instead of eq.(JT9l). To satisfy this condition we assume the following relations between £( m ) and 

^ m \\,fi;t) = r ] ( m \\,v)*^ 1 \v^;t), m > 1, r/«(A, u) = X 2 (A, //), (37) 

where r]^ m ' (A,//) are some n x n matrix functions, which will be specified below. We also 
define t-dependence of £^ as follows: 

dl } (A, n; t) — TH ( A, „) * {v, m *) , (38) 

where T^ m )(A, i>) are some n x n matrix functions. Substituting eqs. (j3"T]) and fl38|) into eq. (l36|) 
we obtain the following representation for a^ m \ m > 1: 

a M(A, /i) = a«(A, i/) * T( m )(z/, z>) * (T^y\u, ft) * faW)" 1 ^, /x). (39) 

In turn, the compatibility condition of eqs. fl38|) requires 

T (m) ^ T (n) _ T (n) ^ r (m) = 0, (40) 

T (m) (A,/i) = ^T(A, z/)r (m) (z/)j ^T _1 (z/,/i), [r( m )(z/), r^(z/)] = 0, V n, m, (41) 

where T(A, /i) and T^(fi) are some no x no matrix functions of arguments. Thus compatibility 
condition of the system fl34l generates eqs. ()37H4Tj) . 
Now, integrating eq. f|38|) with m = 1 we obtain: 

e (1) (A, u; t) = (T(A, ^)e^ * W M**) * T -1 (i/, z>) * £ (z>, «)■ (42) 

Finally, integrating eq. f)34p with m = 1 we derive the following explicite formula for ^(A, H]t): 

X (A, u; t) = a« (A, z/) * (T«)~>, z>) * u; t) + Xo (A, «), (43) 

where Xo(A,/i) is n x n matrix integration constant. 

It is convenient to rewrite expressions for and i (m) (see eqs. (j3"Tl) ) using eq. (l39l) as 
follows: 

A (m) (A,u) =a (1) (A,z/)*r (m) (^,"), i (m) (A) =a (1) (A,z/)^f (m) (^), (44) 
r( m )(A, ^) = T( m >(A, i/) * (tW)" 1 ^, i/) * faM)- 1 ^, A) * /3 (m) (A, u), 
fM(A) = T (m) (A, i/) * (T^)" 1 ^ i?) * (f/" ) -1 ^, A) * ^ (m) (/i)- 

Since r/ (1) (A,u) = X 2 (A,u), one has = and f (1) = Representations (SI) will be 

used in the rest of Sec J2.31 
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2.3.1 Internal constraints for r( m ), T^, T and r< m > 

Remember that definition of ^ m ' ) in terms of x> see ec L- (|35|) . must be consistent with eg. ( 143]) . 
This requirement generates a set of constraints on T^ m \ T^ m \ T and r^ m \ To derive these 
constraints, we apply operator (fi( m > + 0^ m 'P)* to eq. p3|) from the left. One gets 

f M _ = (p(m) + ^Np) # a (l) ^ ( T (l))-1 * £(1) + (£(m) + ^Hpj # ^ (45) 

Substituting eq. (!37|) for £( m ) one has to obtain an identity valid for any This requirement 
yields, first of all, the following equation relating (3^ m \ and Xo (the first constraint on the 
functions r( m ) and f ( m >): 

pim) + p(m)p^ ^ ^ + ^(m) = Q ^ (46) 

r< m >(A, v) + f ^(X)P(u)) * xo(^,/i) + T^(u^) = 0, (47) 



and the following definition of -q^: 

r] (m) = [p^ + /3 (m) p) * * (T^)- 1 , m = l,...,P>. (48) 

Finally, eq. f|48|) written in terms of and r( m ) yields: 

T( m )(A,/i) = T«(A,z/)* (r( m V,^) + f( m V)P(i>)) * * (Tt 1 ^ 1 ^,//), (49) 

m — 1, . . . , D, 

where are represented by eq. fT4T]) in terms of T and r^ m \ This is the second constraint 
imposed on the matrix functions T^ m \ T^ m \ T and r^ m \ 

Constraints obtained in this section are produced, generally speaking, by system (fl5|) and 
its compatibility condition. For this reasong we refer to them as the internal constraints. In 
contrust, the external constraints will be introduced "by hand" in order to derive the nonlinear 
PDEs, see Theorems 2.2 and 2.3. 

2.3.2 First order nonlinear PDEs for the functions V^ m '(t), m = 1, . . . ,D 

Theorem 2.2. In addition to eqs. (Ill P15j) . we impose the following external constraint for the 
matrix functions A( m )(A,yu) and A^^X): 

D D 

^4 (ra) (A ) i/)*i ( " ) (i/)B w = ^i (m) (A)P ( ""' ) , n=l,...,D, (50) 

m=l m=l 

where p( mra ) and p( mra ) are some n x n constant matrices. Then no X no matrix functions 
V^ m \t), m — 1, . . . , D, are solutions to the following system of nonlinear PDEs: 

D 

Yl [(^ + V {m) V {n A B (mn) + V (m) ^A {n) B {mn) + P (mn) ) =0, n=l,...,D, (51) 

m=l 

A&> = P(A) *i (n) (A). 
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Proof: Applying operator *A( n > to the eq.( fT6l) from the right one gets the following equa- 
tion: 

E (mn) ^ ._ y(n) ^ + y( m ) ^ fy(n) ^ + ^(n) ^ + jjimn) ( t ) = Q, (52) 

which introduces a new set of fields U^ mn \ 

U imn) (t) = (w(/i;t) - P(jij\ *A' m »(M,A)*I (T!) (A). (53) 

Due to the relation (|50|) . we may eliminate these fields using a proper combinations of eqs. fl52l) . 
Namely, combination Y^L=i E ( - mn \t)B ( - mn '> results in the system ( l5"B . ■ 

Reduction 1. Let 

p(mn) _ _j\riQ{mn) ^ 

Then eq. tloTj) reduces as follows: 

D 

+ V (m V (n) ) B (mn) = (55) 

m=l 

Note that eq. (!55|) coinsides with the linearizable eq.( )29|) if B^ mn ' = B^. However, eq.( )29|) is 
supplemented by constraints (I25p . which are not valid for eq. (1551) in general. Constraint ( I50p 
reads in this case: 

D 

^{A {m \\v)+A^ m \\)P{v))*A {n \v)B { - mn) ={), n = l,...,D. (56) 

m=l 

Reduction 2. Let 

f (m) (A) = f(A)C (m) , C (1) =/ no , C^bW = B {m \ m,n = l,...,D, (57) 

in addition to reduction ( I54p . Here C^ 71 ' and B^ are some Uq x no constant matrices and T(A) 
is no x n matrix function. As a consequence, we obtain 

IM(A) = i(A)C (m) , i(A) = a (1) (A, u) *f(u), V (m) (^) = V(t)C (m) , f(A) = f«(A). (58) 

Then the system of .D equations f )55|) reduces to the following single PDE: 

(V^ + VC^V) £( m > = 0. ( 59 ) 

m=l 

This equation is written in the Introduction, see eq.([T]), and may be referred to as a multidi- 
mensional generalization of (l+l)-dimensional C-integrable iV-wave equation (|30|) . Reduction 
fl5T|) changes internal constraints ( 14"T|) and (l4"9~j) as follows: 



Y^(A, i/) + f(A)C< m )p(z/)) * X o{u, u) + r^(A, u) = 0, (60) 

T(™)(A, fi) = TW(A, i/) * (t^\u, z>) + f C (m V)^)) * a (1) (z>, At) * (T {1) )-Hfi, At), (61) 
m = 1, . . . , D. 

In turn, external constraint ( 1361) reduces to the following single equation: 

D 

^(A (m) (A, v) + i(A)C (m) P(z/)) * i(z/)P (m) = 0. (62) 



m=l 
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2.3.3 Second order nonlinear PDEs for the functions V^ m \ m = 1, . . . , D 

Theorem 2.3. In addition to eqs. (lll|ll5l) . we impose the following external constraint for the 
matrix functions A( m )(A,/i) and A^ m \\) (instead of constraint (1501) ): 

D 

A {m) (A, v) * A {n) (u, n) * A® (n)B {mnl) = (63) 

m,n=l 

D D 

J2A {m) {X,n)*A {n \fi)P {rnnl) + ^i (m) (A)P (m °, 



m=l 



where B^ mnl \ pi™*), and P (mZ) are some constant no x no matrices. Then no x no matrix 
functions (t) are solutions to the following system of nonlinear PDEs: 

D 

Y WS> m + iy {n) y {l) )t m + v {m) v^ + v {m) v {n) v (v> B {mni) + (64) 

m,n=l 

D D 

{V^ + V^V^A^ + Yl y {m) ^t l) = 0, l = l,...,D, 

m,n=l m=l 

D 

j[mnl) _ g(mnl) _|_ p(mnl) ^( m = J[(n) pjmnl) _ p{ml) 

n=l 

Proof: Applying operator to the eq. (|T6l) from the right one gets eq. (1521) . which 

introduces a new set of fields U^ mn ' > defined by eq.( l53l) . However, we may not eliminate these 
fields from the system fl52l) because constraint (1501) is not valid in this theorem. Instead of 
this, we derive nonlinear PDEs for these fields as follows. Applying operator *A^ n ) * A^ to the 
eq. ffTB"]) from the right one gets 



E {mnl \t) := U%!\t) + V {m \t)U {rd) (t) + U {mnl \t) = 0, (65) 

where one more set of matrix fields appears: 

U {mnl \t) = (W(A;t) -P(A)) * A (m) (\,u) * A (n) (i/,/i) * A (l \fx). (66) 

Due to the constraint (163]) . these fields may be eliminated in a proper combination of equations 
(165|) . namely, Emn=i^' mnl '^' mn ^ Then, substituting U ( - mn '> from eq. (l52|) one ends up with 
eq.@3). ■ 

Emphasize that nonlinear equations (1511) and (164"1) do not represent commuting flows since 
we assume different external constraints (1501) and (1631) for matrix functions A^(A, ji) and 
A( m )(A) deriving these equations. These two constraints are not compatible in general. 

Reduction 1. Let 

p{mnl) _ p{mnl) ^ 

D 

p(ml) _ 'y ^ j^(n) p(mnl) 



n=l 
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Then eg. AM]) reads 



D 

[ V tlL + (V in) V (l) )t m + V™vg + V {m) V {n) V {l) B {mnl) = 0, I = 1, . . . , D. (68 



m,n=l 

Constraint fl63|) reduces to the following one: 

D 



( y A {m X\,v)*A {n \is, f i) + A {m \\,is)*A( n \is)P( f i)+ (69) 

m,n=l 

i (m) (A)P(z/) * i (n) (z/)P(^)) * A il \n)B {mnl) = 0. 



Reduction 2. Along with reduction (J6"T|) we consider reduction (I57II58P with 

£<(0 g(mnl) _ g(mn) ^ 

where p( mn ) are some constant matrices. System (168]) reduces to the following single PDE: 

D 

i Vt ^ + ( Vrc(n)y )^ + VC im) V tn + VC (m) VC (n V] P (mn) = 0. (71) 



m,n=l 



This equation is written in the Introduction, see eq.(j9]). Internal constraints f l60|) and (ISTj) 
remain valid for this well. External constraint ( 1691) reduces to the following single 

equation: 

^ ^A {m \X,u)*A {n \u,fi) + A im \X,u)*A(u)C {n) P{fi)+ (72) 

m,n=l 

i(A)C (m) P(z/) * A{u)C {n) P{fi)) * A{fi)B {mn) = 0. 



2.3.4 Solutions to the first order nonlinear PDE ( 1591) 

The problem of richness of the available solution space will be considered for the nonlinear PDE 
(R)9~|) . We show that solution space to this equation may be full provided that all constraints 
(16"U|) . (IBTj) and fl6"2l may be resolved keepeng proper arbitrariness of functions r^ m \v). Examples 
of particular solutions will be considered as well. 

On the dimensionality of the available solution space. We estimate the dimensionality 
of solution space for small \- I* 1 this case formula (TTTT) yields W(X; t) ps P{y) * x{ v -> A; t) an d 
formula (TTTjl gives us 

V(t) « (P(z/) * x(^, A; t) - P(A)) * A(X). (73) 

By construction, if all r M(!/) (m = 1,... , .D) are arbitrary functions, this expression preserves 
the following arbitrary n X n matrix function of all D variables: 

F(t) = P*a (1) *(T (1) )" 1 *£ (1) *i = J g(is)e^ T( ' Hu)t 'g 2 (is)dn 2 (is), (74) 
^(i/) = P{v) * a (1) (^,/i) * (T^)- 1 ^, A) *T(A, i/), ^ 2 (i/) = T-^i/, A*) * 
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However, dimensionality of this function reduces due to the presence of constraints ( 1601) . ( l6Tj) 
and (162]) which impose relations among T^ m \v). An important question is whether the dimen- 
sionality of the function (1711) may be equal to D — 1, which is necessary for fullness of the 
solution space. At first glance, we may expect the positive answer. In fact, eq. ( 160]) may be 
satisfied using special structures of T^, P and Xo-, as it is done in the example below, see 
eqs.( )89p . Next, eq. (l6ip relates r^ m ' with and r^ m ^ which, in general, keeps arbitrariness 
of all and consequently does not restrict dimensionality of the above written arbitrary 

function. Finally, eq.( l62p must be considered as a single relation among r^ m \u), m = 1, . . . , 
redusing the dimensionality of the function (|74p from D to D — 1, which means the full di- 
mensionality of the solution space. Thus, we may expect new completely integrable nonlinear 
PDEs in the derived class of equations. 

We have outlined a rough analysis of the solution space dimensionality. More detailed 
analysis must be carried out for particular equations and remains beyond the scope of this 
paper. 

Construction of explicite solutions. Now we derive a family of particular solutions in the 
form of rational functions of exponents. Solitons and kinks are the most famous representatives 
of this family. To derive such solutions, we take 

M N 

dfii(A) = ^ 5(A - a,i)d\, d0 2 (A) = J2 S ( X ~ 6 *) dA > T ( x > =MK v), (75) 

i=l i=l 

2^1 — > Imu , %2 Inu - 
Then all integral equations reduce to the algebraic ones. We use notations 

X(ai,ai) ••• x(ai,a M ) 



W = [ W(ai) ■ ■ ■ W(a M ) ] , X = 
Xo(ai,ai) •■■ Xo(ai,OAf) 

Xo{ a M, Oi) ■■• Xo(. a M,0>M) 



Xo 



x( a M, ai) • • ■ x( a M, o>m) 



(76) 



Of. 



(i) 



a {1) {a M ,bi) 



a (1) (ai,67v) 
a w (a M ,b N ) 





" rW(6i,oi) ■ 


■ rW(6i, flJf ) " 




" f(6i) " 








, f = 






_r^(b N , ai ) ■ 


• r( m )(6^,a M )_ 




. f (b N ) _ 





£o(&i,ai] 








_ 6o(&iV, a l 


f m ) = 


diag(r (m) 



Co(&l,«Af) 



P=[P(a 1 ) ••• P(a M )], 



Solution V is given by eq.( lT7|) together with reduction ( l57|58l) as follows: 

y = (w-P)A 

where W is solution to eq. ffTTj) : 

W = Px(X + hln )- 1 



(77) 



(78) 
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Substituting eq.liTSl) into eq.(l77l) we obtain 

V = P(x(x + W 1 - /Mno)a {1) r. (79) 
Since \ is given by eq. (H3l) . one has 

X = d«(f«)- 1 e ^^6 + Xo, (80) 

where we substitute eq.( H2|) for ^\ 

Matrices T^ m \ r and f*" 1 ) must satisfy constraints ( 16"Uj) . flBTj) and (I6"2"j) . which read in our 
case, m = 1, . . . , D: 

(T (m) + f C (m) P) xo + f M = 0, (81) 

f(m) = f (l) ^f(m) + f C (m)p^ £(1)^(1))-^ (g2) 

^ a (1) (f (m) + rC (m) P)a (1) f\B (m) = 0. (83) 

m=l 

Analisis of eqs.f l8Tj) points on two different types of solutions to them. First type is assotiated 
with det(xo + Iun ) 0- Then eqs. f[8"Tj) may be solved for and one can show that 
multidimensional PDE (|59|) may be splitted into a set of independent compatible Ordinary 
Differential Equations (ODEs). We will not consider this case. Second type is assotiated with 
det(xo + hin ) = and leads to truly multidimensional solutions to eq. (!59|) . Namely this case 
is considered hereafter. 

Looking for the particular solutions to eq. (1811) we decompose it into two equations: 

PXo = 0, f( m )(xo + W=0, (84) 

which means that the rows of P and f ^ are orthogonal to the columns of Xo an d Xo + Imu 
respectively. 

Note that eq. fl82l) with m = 1 reads: 

f (D = (r«+fp)a«. (85) 

Since, must be invertable, we require M > N. Then eq. ([8"3"j) may be simplified removing 
as a left factor in this equation: 

D 

J^(fM + f C (m) P)d (1) f P (m) = 0. (86) 

m=l 

Simple example of solution. We consider the three-dimensional nonlinear PDE fl59|) . i.e. 
D = 3: 

3 

( Vt ™ + VC [m) V) B{m) = 0- (87) 

m=l 
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Let n = 2, M = 5, N = 2, 

B (i) = C {i) = j 2j B (i) = diag ( 6 « &«), C« = diag(cf ) , cf), i = 2, 3. (88) 

In order to satisfy eqs. (!82|) . (|84|) and (!86|) we take the following matrices Xo> A f f : 



Xo 



(i) 



— -^2 ^2,8 

-^2,2 -^2,8 

' h 

h 

Jo 



F 2j 2 — Z, F 2j 8 — [J a Z J5 Z], 

io = (« (1) ) T 



P = [Z Z Z J x Z] , 

J 2 Z Z Z Z 

z z z z z 



f« = 
f = f w 



z 

■h 



r {2) = f w = z 4 , 10 , 

f (1) = diag(l, 2, 3, 4), f (2) = diag(0, 0, 3cj, 2) , 4cS 2) ) 



f< 3) = diag(0,0,34 3) ,4ci 3) ). 
In addition, we obtain expressions for c\ , i = 1, 2: 



(3) 



l + 6i a) cj a) 



z = 1,2. 



Here Zjj and Z are z x j and 2x2 zero matrices respectively, 
J = [/ 2 Z], J x 



' 


1 " 


, J2 = 


' 1 


" 




" 


3 " 


1 








2 


4 






J, 



8p 3 + £2(6 - p 3 ) - 4p 4 - pi(12 - p 4 ) 

2(6pi - 3p 2 - 2p 3 +p 4 ) 
(-12 + 3p 2 -p 4 )(-2p 3 +p 4 ) 

12(6pi - 3p 2 - 2p 3 + p 4 ) 

x 4(3pi - p 3 ~ 6) " 

12 - 3p 2 + p 4 





1 




, , P2 _ n 
4 12 



(89) 



(90) 



(91) 



(92) 



Diagonal elements of matrices B«, i = 2,3, and C^ 2 -* remain arbitrary. Substituting eqs. (!89p 
into eq.( l79ll80l one obtains V as a rational expression of exponents: 



V(t) = I 



-4(e> 2 +p 4 ) /i.e" 1 
/ 21 e" 2 -3(e> 3 +P4) 
£) = e 771 ^! + e m p 2 + e^pa + p 4 , 

, = 3 (-12 + 3p 2 - p 4 )(p 2 p 3 ~ P1P4) , = 16 jjgi ~ 6 P2 ~ 4p 3 + 2p 4 
112 4(-12 Pl + 6p 2 + 4p 3 - 2p 4 ) ' /21_ ' -12 + 3p 2 -p 4 



(93) 
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Here 

77! = 4ti + 4ci 2) t 2 - 4t(! + &iM 2) )*3, (94) 

m = 3tx + 34 2) t 2 - + 4 2) 4 2) )t 3 . 

We see that all elements of V are kinks. 

Relations (1901) show that not all matrix coefficients in eq. fl87p are arbitrary. They are related 
by the following equation: 

C®B® + C®B®+I 2 = 0. (95) 

Thus we have constructed particular solution to the three-dimensional nonlinear PDE (1BT1) with 
diagonal matrices related by eq. fl95|) . 

3 Multidimensional generalization of -S-integrable PDEs 

3.1 Starting equations 

Algorithm developed in this section is based on the same equation (iTTj) with different function 
x(t), which is defined by the following system of equations 

Xt m {\^t) = (A^\X,u)+A^\X)P(u))*x(^^t)-x(X,u;t)*A^\v, t i), (96) 
m — 1, . . . , D 

instead of system (fl5|) . Here, again, A^ m \X, v) and A^ m \X) are n$ x no matrix functions of 
arguments. 

Matrices A^ and A^ have to provide compatibility of system (196]) . Similar to eq. (fT5l) . 
there are two different methods that provide this compatibility. The first one yields the classical 
S-integrable nonlinear PDEs, Sec J3.2} while the second method yields a new type of nonlin- 
ear PDEs whose complete integrability is not clarified yet, Sec J3.3l However, our algorithm 
supplyes, at least, a rich manifold of particular solutions to these PDEs. 

3.2 First method: classical ^-integrable (2+l)-dimensional iV-wave 
equation 

Consider the compatibility condition of eqs. fl96l) in the following form: 

[A^ (A, y) + i (m) ( A)P(i/)) * X t n (v, K t) - xu (A, v\ t) * A^ (u, = (97) 
(AW(A, u) + I w (A)P(i/)) * xt m {v, K t) - XtjX, v- 1) * A™(y, fi), V n, m 

Substituting eq. fl96l) for derivatives of x we reduce eq. fl97|) to the following one: 

(L im) * L (n) - L (n) * L (m) ) * x - X * {A [m) * A {n) - A {n) * A {m) ) = 0, (98) 
L im \X^) = A^ m \X,/i) + i (m) (A)P(/i). 
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Since eq.( l98|) must be valid for any function xif) (which is a solution to the system ()96l)). it is 
equivalent to two following equations relating matrix functions A^ m \ and P: 

A {m) (A, v) * A (n) (u, n) - A (n) (A, v) * A (m) {y, n) = 0, (99) 

^ra) # jjn) _ L (n) # L (m) = q 

L (m) (A, z/) * i (n) (z/)P(^) - L (n) (A, i/) * A {m) (u)P{>i) = (100) 
i (n) (A)P(z/) * A (m) (z/,/i) -i (m) (A)P(z/) * A (n) (z/, /i). 

In order to satisfy eq. fllOOp we require the following representation of A' 771 ' (A): 

A (m) (A) = A(A)P (m) , [P (m) ,P (n) ] = 0, (101) 

where A(\) and p(' m ) are no x uq matrix function and constant matrix respectively. Then 
eq. (11001) is equivalent to the following system: 

L (m) (A, v) * A(z/)P (n) - L (n) (A, v) * A(z/)P (m) = 0, (102) 
P (n) P(z/) * A {m) (u, /i) - B {m) P(u) * A {n) (y, /i) = 0. (103) 

Eqs. fl99lll02lll03p represent three constraints for matrices A^ and A. 

Theorem 3.1. Let the matrix function W(X;t) be obtained as a solution to the integral 
equation (TTTj) with x defined by eqs. (!96|) supplemented with constraints (I99II102|I103I) . Then 

1. Function W(X; t) satisfies the following system of compatible linear equations 

P> m >(M := B^(W t Jfi;t) + V(t)B^W(^t) + W(u;t)*A^\u,ii))^ (104) 
B^ (W u {^ t) + V(t)B^W(ii; t) + W{v- 1) * A^{u, n)) = 0, 
V(t) = (W(u ] t)-P(u))*A(u), (105) 
n, m = 1, . . . , D. 

2. Matrix field V(t), given by eq. (I105j) . satisfies the following S'-integrable A^-wave equation: 

( B(n) ( V tm + VB {m) V) - B^ m \V tn + VB^V))bV = 0, (106) 

perm(n,m,l) 

where perm(n, m, I) means clockwise circle permutations. 

Proof: 1. To derive eq. (1104[) . we differentiate eq. ffTTj) with respect to t m . Then, in view of 
eq. f[9"6"j) . one gets the following equation: 

( M ; t) := P(A) * A< m > (A, v) * x (y, fi; t) = P» (i/; t) * A*) + M", /*)) , (107) 
(w t) = W tm (//; t) + (t)W(K t) + W{y- 1) * A< m > {v, n) . 

Due to the constraint (11031) . LHS in eqs. fll07p may be removed using the following combination 
of these equations: 

B (n) S (m) _ B {m) £ (n) ^ (108) 

(P (n) P (m V;t) - B^E^ n \u;t)) * (x{v,Wt) +Zi{v,l*)) = 0. (109) 
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Since operator *{x(v, A*) + Xi(z/, is invertable, one gets 

£(nm) . = B (n) E (m) _ fl N £ («) = Qj (HO) 

which coinsides with eq.l l 104ft . 

2. In order to derive nonlinear PDE (11061) we consider the following combination of eqs. (11041) : 

E (nm) *AB {l) , (111) 

perm(n,m,l) 

which yields eq. (I106p in view of constraint (11021) . ■ 

S'-integrable PDEs will not be considered in this paper. 

3.3 Second method: new class of nonlinear PDEs 

We will use double indexes hereafter in this section, i.e. 

D D 2 Dx 

y(m) = v (m ima)j tm =t mim2} J2f {m) =H Y,f {mim%) > V /> D = (D 1 ,D 2 ), (112) 

m=l rri2=lmi=l 

and notation p 1 ' = f^ u \ V/, for the sake of brevity. For instance, = r^ n \ r/ 1 ) = r]^ 11 ' and 
so on. Similar to Sec J2.3[ we will use representations fl3T|) for the matrix functions A^(X, /i) 
and A*- m ^(A) and write eq.( |96l) in the following form: 

Xt m (\v,t) = a (m) (A,^)*e (m) (^/i)-e (m) (A^)*/3 (m) (^^), (H3) 
£ (m) (A,/i)= ^ m \\,v) + fr m \\)P{v)) * X (v,r,t), (114) 

e (m) (A,/i) = x(A^)*a( m )(z/,/i), m i = l,...,D i , z = l,2. (115) 

Then the compatibility condition for system (196]) is equivalent to the compatibility condition 
for system (11131) which reads 

a (m) ^ £m) _ gm) ^ ^ = a (n) ^ ^(n) _ gn) ^ ^ y ^ m _ (ng) 

To satisfy this condition we, first, assume that ^ m - ) and <5^ m ) are expressed in terms of £^ and 
£W as follows: 

^ m \X, f i ] t) = ^ m \X,v)*^ 1 \u, f i), m 1 + m 2 >2, r/^A, z/) = Z 2 (A, z/), (117) 
^ m \\, f i;t) = ^ 1 \\,v;t)*r ] (m \v, f i), m x + m 2 > 2, ijWfo //) = X 2 (A,i/), (118) 

where r]^ m \u, /x) and fj^ m \u, jj) are some n x n matrix functions which will be specified below. 
Second, we define £^ and £f in terms of £^ and as follows: 

^(A,/i;t) = T^A.^^W;*), (119) 
^(A,/*;t) = PHA,M)*T( m W), (120) 

where T( m )(A, v) and T( m )(A, z/) are n x n matrix functions. At last, substitute eqs. (lll7H120|) 
into eq. fll 16|) . Since resulting equation must be valid for any possible £^ and we get the 
following expressions for a (m) and /3 (m) : 

a (m) (A,/i) = (a (1) (A,z/)*T( m )(z/,i>)*(T( 1 ))- 1 (z>,z/)) * (t/^)" 1 ^), (121) 
/^(A,^^)- 1 ^^ (122) 
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In turn, compatibility of eqs. (I119j) and f II 2 f) requires 

T M^ T (n)_ T (n)^ T (m) =Q ^ T (m) (X, fi) = (T (X, u)t^ 1 \u)) * T' 1 (u, fi) , (123) 

fH^(»)_f(n)^H =0 f (m) (A,/i) =f- 1 (A,z/)*(f (m) (z/)T- 1 (z/, y u)), 
[ r M (zy ) ) r (») (!/)] = 0, [f < m) (i/) , (i/)] = 0, 

where T(A, z/), T(A, z/), r^ m ^(z/) and f^(zv) are no xno matrix functions (compare with eq.( l4Tj) ). 
Thus, compatibility condition of eqs. (H13[) generates eqs. (II 17111231) . 

Now we may integrate eqs. (|119|) and (I120p obtaining the following expressions for £ < - 1 - ) and 

(124) 

£ (1) (A )M ;t) = (T(X,u)e^ (i) ^) ^T" 1 ^, v) *£ (z>,/i), , (125) 
= £o(A,z,)*T->,z>)* ( e E£*' w (**r (z>;/i) ) . (126) 



Here £o(^)/f) an d £o(A, are % x matrix functions. Finally, integrating eq. (H13|) with 
mi = Tn2 = 1, one obtains 

X (A, /i; f) = a« (A, i/) * (T«)->, v) * £ (1) (*>, " (A, * (T«)->, *) * /3 W A*)- (127) 

It is convenient to rewrite eqs. (l3Tj) using eqs. (1121ft and (1122ft as follows: 

A (m) (A,/i) = a^\\u)*^ m \u,u)*fi x Xu, t i), (128) 
IM(A) = a (1) (A,z/)*f (m) (z/), Vm, (129) 

where 

7 M (A,aO = T^\\v)*{T^y\vrv){ri m \\v^)* (130) 
(^))- 1 (/i, / i)^(f«)- 1 ( / i,A)^f( m )(A, / i), 

f(m) (A) = T (m) (A,Z/)*(T (1) )~ 1 (z/,Z>)*(77 (m) )- 1 (z>,/i)*^ m) (/i). (131) 

Representations (j!28|129p will be used in the rest of Sec l3.3l 

3.3.1 Internal constraints for f( m \ T, f, r( m ) and r( m ) 

Note that the definitions of and in terms of Xi i- e - eqs. ( 11141) and ( j 1 15j) . must be 
consistent with eq. ( 1127ft . This requirement generates some constraints for a (1) , f^ m \ T, 
T, and f( m \ To derive constraints assotiated with definition ( 1114ft . we apply operator 
(/?( m ) -f p^P)* to the eq. (1127ft from the left obtaining the following equation: 

£ (m) (A,M = (/3 (m HA^) + ^ (m HA)W)*^ 

(^ m \X, v) + fr m \X)P{v)) * £«(i^) * (rW)" 1 ^, p) ★ ai). 

Substitute eq. (11171) for and require that resulting equation is identity for any Then 
eq. (1132ft becomes decomposed into two following constraints: 

(/3( m )(A, !/) + ^(A)P(z/)) * £ Q (y, y) = 0, (133) 
V {m \X,fi) = (/3 (m) (A, i/) + /3 (m) (A)P(i/)) * a (1) (i/, z?) ★ (T^)-\u, //), (134) 
mj = l,...,A, i = l,2. 
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Similarly, to derive constraints assotiated with definition ( 11151) . we apply operator to 
the eg. (11271) from the right. One obtains 

£M = a « * (TW)- 1 * «M - * (T^)" 1 * /3« * a H. (135) 

Substitute eq. ffTT8|) for and require that resulting equation is identity for any Then 
eq. (11351) becomes decomposed into two following constraints: 

e (A,z/)*a«( I /,/i) = 0, (136) 
r/ (m) (A,/i) = -(fW)- 1 (A,i')^ (1) (f,i')*a w (5 ) /'). (137) 

Eqs. (ll33p and (11361) represent two constraints for a^ 1 ', f3^ m ' ) and (3( m \ Eqs. (ll34p and (11371) with 
mi + m 2 > 2 may be considered as definitions of r]^ m ' and fi^ m \ However, since r)^(\,fj,) = 
fj^(X, /i) =X 2 (A,/i), eq. (11371) with mi = m 2 = 1 yields: 

f( 1 )(A,^) = -/3«(A,z/)*a«(z,, / i), (138) 

In turn, eq. (H34p with m\ = m 2 = 1 in view of eq. (H38p yields 

T( 1 )(A,/i) = -f( 1 )(A, y u) + ^ 1 )(A)P(z/)* a ( 1 )(z/, / i), (139) 

Both eqs. (U38p and (I139P may be treated as constraints for and 0^'. 

Now we simplify eq. (11371) . mi + m 2 > 2, replacing (3^ * in accordance with eq. (11381) . 
One gets 

fjW+^m) = T (m) + (T (l) r i e q-pH]) (14Q) 

7 (m) (A,^) = (f^)-\\,v)*f( m \v^), mi + m 2 > 2, (141) 

which is the definition of 7^ . Deriving eq. ( 11411) we assume invertibility of the operator -kT^ m > , 
V m. 

Constraints (1 133j) and (11341) may be simplified multiplying them by fj^ from the left, using 
definitions of (jI2B . (3^ fTl22l . f ^ (11311 and eliminating with eg- flT^O]) . We 

obtain in result: 

(f «)- 1 (A, 1/) *f <»>(„, z>) *0W(u,(i) + f^(X)P(fi)) * fo(/i,M) = 0, V m, (142) 



^(m) * JWp + a (i) = T M + T H (143) 
f (m) (A)P(z/) * a (1) (z/, /i) = T (m) (A,/i) + T (m) (A,/i), m x + m 2 > 2. (144) 

One has to take into account that T^ m \ T( m ' are represented by eqs. (11231) in terms of T(A,/z), 
T(A, /i), r( m )(A, /i) and f ( m )(A, /i). Thus, we have obtained a set of constraints for the functions 
aW, f ( m ), T, f , r( m ) and f^: eqs. (ll36m380139m42m44l) . 

Note that, similar to Sec J2.3.1[ all these constraints are generated by system (1961) and its 
compatibility condition. For this reasong we refer to them as the internal constraints in order 
to defer them from so-called external constraints wich will be introduced "by hand" for the 
purpose of derivation of the nonlinear PDEs, see Theorems 3.3, 3,4, 3,5. 
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3.3.2 System of compatible linear equations for W(X;t). 

Theorem 3.3. Let matrices A^ m '(X, /i) satisfy the following external constraint: 

Di 



L^P(X) * A (m) (A,/i) = S {m2) P(X), (145) 



mi=l 



where L^' 171 ^ and S^ 7 ™ 2 * 1 are some n x n constant matrices. Then matrix function W(X;t) 
obtained as a solution to the integral equation (fTTj) with x defined by eq. ([9"61) is a solution to 
the following system of compatible linear equations 

E (m2 \X;t) := (146) 

Di 

L {mi) {W tm (X; t) + V (m \t)W(X; t) + W(ji; t) * A^(ji, X; t)) = S^W^X; t), 

mi=l 

mi = 1, . . ., D h i = 1,2 

where is given by eq. f|T7|) . 

Proof: To derive eq. f)146p . we differentiate eq.f llip with respect to t m . Then, in view of 
eq.(J2SJ, one gets the following integral equation: 

£ {m \fi;t) := P(u) * AW(y, X) * x(A,//;t) = E^\u;t) * { X {u,^t)+X x {u,^)), (147) 
( M ; t) = W tm (n;t) + {t)W{^ t) + W(u;t)* A™ {u,n). 

Consider the following combination of eqs. using constraint ( 1145 j) , 

one gets in result: 

Di Di 

J2 L ^) £ ™ ;= ^ L^M^f) * (*(!/, A*! *) + Z 1 (i/,/i)) = (148) 

mi=l mi=l 

Since operator t) +Xi(z/, //)) is invertable, eq. (j!48|) is equivalent to eq. f ll46[) . ■ 

Remark: Similar to eq.f lTB"]) . eq. (j!46j) is, strictly speaking, nonlinear equation for W(X;t) 
since V^ m \t) are defined by eq.flU} in terms of W(X; t). 

System (j!46p is an analogy of the overdetermined system of linear equations in the classical 
inverse spectral transform method. According to this method, nonlinear PDEs for potentials 
of overdetermined linear system appear as compatibility conditions for this system. However, 
nonlinear PDEs may not be obtained by this method in our case because of the last term in the 
LHS of eq. (I146j) . Instead of this, we represent another algorithm of derivation of the nonlinear 
PDEs in Secs J3X3l and ETTOl 

3.3.3 First order nonlinear PDEs for the fields V {m \t), rrii = 1, . . . , A, i = 1, 2. 

Theorem 3.4. In addition to eqs.( lll|l9"6"j) and external constraint (I145j) . we impose one more 
external constraint: 

D 2 D 

A {m \X,u) *i (,l V)fl W = ^A {j \X)P {mijn \ m = l,...,Di, i = l,2, (149) 

m 2 =l j=l 
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where p( m2n ) and p( mi i n ) are some uq x uq constant matrices. Then no x no matrix functions 
V^ m \t) are solutions to the following system of nonlinear PDEs: 

D Di D 

£ j ( m i) -f y( rn )y( n ) _|_ yM^fa) J flfrnan) _|_ £( mi ) yCOpOij™) _ (150) 

m=l mi=l J=l 

£>2 

^ S h) y W i j(m 2 n) i = P(A,z/)*i(")(z/), n 4 = l,...,A, i = 1,2. 

7713=1 

Proof: Applying operator *y4^ n ^ to eq. (11461) from the right one gets the following equation 

E (m 2 n)^ = E (m*)( X . q * ^ n \\) := (151) 
mti m t2 

7711=1 7712=1 

which introduces a new set of fields U^ mn \ see eq. fl53|) . Due to the constraint (I149p . we 
may eliminate these fields using a proper combinations of eqs. (11511) . Namely, combinations 
^^ i£ (m 2 «) fi (m 2 n) results in the system ([150]). ■ 

Reduction 1 The derived equation ( 11501) admits the following reduction for its coefficients: 

pirrnjn) = -A^R^5 jimi , S {m2) = 0. (152) 

Then eq. (11501) reads: 

D 

L{mi) + V {m) V {n) ^j P (man) = 0, (153) 

m=l 

which is represented in the introduction, see eq.©. This reduction does not effect constraint 
(11451) while constraint (11491) reads: 

d 2 

^(A (m) (A,z/) + i (m) (A)P(zy))*iW(z/) J R (m2n) =0, n i = l,...,A, i = l,2. (154) 

7712=1 

Reduction 2 Along with reduction ( 11521) . we consider reduction ( |57ll58l) with 

C (n) R (m 2 n) = R (m 2 ) _ ^55) 

where and R^ 2 ^ are some no x no constant matrices. Then the system (11531) reduces to 
the single PDE: 

D 

Limi) {Vt m + VC {m) V) P (ma) = 0, (156) 

771=1 

which is written in the Introduction, see eq.©. Reduction f 1 1 5 5 j) does not effect constraint 
(11451) . while constraint (I154p reduces to the following single equation 

d 2 

{A {m) (A, v) + i(A)C (m) P(zv)) * A{v)R (m2) = 0. (157) 

7712 = 1 
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If, in addition, C imim2 \ L (mi) and i? (m2) are diagonal, 

C (mim 2 ) = _ C { m ,) R (m 2 ) L (m 2 ) y = _y D 1 = D 2 = D , (158) 

then nonlinear equation (11561) reduces to the following equation 

V (L {mi) V t L {m2) - L {m2) V t L (mi) + (159) 

m l i m 2 

TTT-2 >)Ti^ 

_^( m i)y(^( m i m 2)y^(m.2) _ £(m 2 )y£i(mim2)^£(mi)\ _ 

which is a multidimensional generalization of the classical (2+l)-dimensional S'-integrable N- 
wave equation (I106p . 

Eq. (ll59p admits reduction 

t mim2 = -ir mim2 , V = -V + , m 2 > m 1 , (160) 

see eq.flH]), which is important for physical applications. 

3.3.4 Second order nonlinear PDEs for V^ m \t), = 1, . . . , D if i = 1, 2 

Theorem 3.5. In addition to eqs. (lllll96p and external constraint (I145P we impose one more 
external constraint on the matrix functions A^ m \X, u) and A^ m \\) (instead of constraint (11491) ): 

D D 2 

A^ m \\,u) *A^(u,fi) *i«(/i) J R( m2n ') = (161) 

n=l mi=\ 

D D 

A (n) (\,fi) * A {p \fi)P (rninpl) + J2A {n) (\)P {minl \ 

n,p=l n=l 

where R( m2nl \ pOi"-p0 ; anc [ p(minl) are some Hq x nQ con stant matrices. Let S^ m ' 2 ^ = for the 
sake of simplicity. Then no x no matrix functions (t) are solutions to the following system 
of nonlinear PDEs: 

D 

L {mi) (u^ l) + V {m) U {nl) + V {m) A {nl) ^ R( m2nl 1 + (162) 

m,n=l 

D 1 / D D \ 

I Jj{np) p{m in pl) _j_ y(n) p(m in l) j _ g^ 

mi = l \n,p=l n=l / 

^)=P(A)*AW(A,/i)*i«(//), /, = 1,...,A, i = l,2, 
where fields [/( mn ) are related with due to eq. (U51j) : 

L{mi) (yt ] + V {m) V {n) + -l/( m )^( n ) + u {mn) J = 0, (163) 

mi=l 

A {n) = P(A) *i (n) (A). 

Proof: First of all, applying operator *A( n > to eq. (11461) from the right one gets equation 
(11631) . which introduces a new set of fields U^ mn \t), see eq. flo^il . This step is equivalent to the 
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first step in derivation of eg. ( 1150jl . However, we may not eliminate these fields from the system 
(I163p because constraint (11491) is not valid in this Theorem. Instead of this, we derive nonlinear 
equations for fields U^ mn \t) applying operator * to eg. (I146j) from the right. One gets 



E (m 2 nl) ^ = ^(m 2 ) ^. Q # A (n) ^ ^ # ftl) ^ ._ 

_^( m i) y[j( nV ) -\- y( m )jj( nl ) _|_ y( m )j[( nl ) _|_ [/( m «0J — Q ; 



(164) 



mi=l 



where the fields ^/( mn are defined by eq. fl66|) . Due to the constraint (11611) . these fields may be 
eliminated taking a proper combination of equations (I164p . namely, X^m 2 =i E^ m2nl ^ R^ m27ll \ In 
result, one obtains eq. (H62p . ■ 



Reduction 1. Let 



p(minpZ) _ r p(n 2 pl) p(mini) _ r p(n 2 l) 



(165) 



p=l 



Then eq. (U62|) reduces to the following one: 



D 

m,n=l 

Zi = l,...,A, 2 = 1,2. 



0, (166) 



Constraint (I145p remains the same (with S^" 1 ^ = 0) while constraint (11611) reads 

D D 2 

J2 Yl ( A{m) ( X ^) + A~( m) (\)P(vf) * (A (n) (z/,/i) * A {l \/i)R {m2nl) - A( n \is)P {m2nl) 



n=l m 2 =l 

mi = 1, . . .,Z>i, /j = l,...,A, 2 = 1,2. 



(167) 



Reduction 2. Along with reduction (I165p we consider reduction (I57II58I) together with the 
following conditions 

jj(ml)(+\ _ Jj( m ) (f^jfji 1 ) ^ (J (l) ft(m 2 nl) _ ^(m 2 n)^ q(u) p(m 2 nl) _ p(m 2 n)_ (168) 

Then the system (11661) reduces to the following single PDE: 

D 

L {mi) (uQ + VC {m) U {n) ^ R {m2n) + (V tm + VC {m) V) P {m2n) = (169) 

m,n=l 



This reduction does not effect constraint (I145p . while constraint (I167p reduces to the following 
single equation: 



D 2 D 



J2 (^ (m) (A, v) + i(A)C (m) P(z/)) * (A (n) (i/, V) * A(v)R [m2n) - A(u)P {m2n ^ . (170) 



m 2 = l n=l 



Emphasize that systems (11501) and (I162|ll63p do not represent commuting flows since con- 
straints (11491) and (I16ip are not compatible in general. 
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3.3.5 Solutions to the first-order equation (11561) 



On the dimensionality of the available solution space. The important question is 
whether the derived nonlinear DPEs are completely integrable. In other words, regarding 
eq. (11561) . is it possible to introduce a proper number of arbitrary functions of P>\D 2 — 1 vari- 
ables in the solution space. 

In order to clarify this problem, first of all, let us consider small x- Then, using equation 
ffTTj) . we approximate V by eq. fl73|) . If t^(u) and -H m )(z/) are arbitrary functions of arguments, 
then the above expression for V involves two arbitrary functions of D\D 2 variables t m , m% = 
l,...,Di,i = 1,2: 

F x (t) = P * * (TW)- 1 * *f = J gil (u)e^ T(l)u g l2 (u)dtt 2 (v), (171) 
F 2 (t) = P * * * (T (1) )- x * * * f = J g 2 i(iy)e^ T( '' >u g 22 (u)dn 2 (u), 

where 

^^P*^ 1 )*^ 1 ))" 1 *^ g 12 = T- l *^*a^*T, (172) 
g 21 = P* * £ * g 22 =f* (T * * * f. 

However, not all 7~( m '(i/) and r' m '(i/) are arbitrary since we have to resolve constraints 
fll36ll38ll39ll42|144ll45lTH7l> . Constraints ffTMT42l may be satisfied using a special structures 
of aS l \ f}( x \ P, £o, £o (see example of explicite solution, eqs. (H84p ). which does not reduce the 
dimensionality of the solution space. Constraints (11381) may be considered as constraints for 
o/ 1 ) and (3^\ Constraints (j!39jl44p relate with and T, so that, generally speaking, 
only one of functions (I17ip remains arbitrary. Two remaining constraints (I145II157P introduce 
D 2 and Di relations among parameters f( m \ = 1, . . . , Di, i — 1, 2, which, in general, reduces 
significantly the dimensionality of the solution space. However, this depends on the particular 
choice of the coefficients and C^ m \ 

This is a preliminary analysis which suggests us to look for examples of completely integrable 
PDEs in the derived class of new nonlinear PDEs. 



Construction of explicite solutions. We construct explicite solutions in the form of ratio- 
nal functions of exponents for the first order nonlinear PDE (I156p with dQi(u), i = 1,2, given 
by eqs. (175|) . Along with notations (176!) we use the following ones: 



£o(ai,&i) 



£o(ai,b N ) 



f (m) =diag(r( m ) (&!),■•• ,r (m) (6jv)), 
T(A, p) = f (A, jS) = X 2 (A, /i). 



(173) 



Field V(t) is represented by eq. ([79l where function x is given by eq. (H27p . In turn, the functions 
£W and £W are given by eqs. (11251) and (11261) . In result, using notations (I76II173I) . we obtain 
the following formula for X- 



(174) 
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Finally, constraints (I136II138|I139II142|I144|I145II157P must be satisfied, which read in view of no- 
tations (I76|173p as follows, m; = 1, . . . , D iy i = 1,2: 



e A (l) = o, 



T 



(1) 



(l)A(l) 



f Pa' 1 ' = f ( - m ' ) + f' m ' 

(f(l))-lf(m)£(l) + fcMpA^ = 0> 



^ L(^)p a (i)(f(i))-if(^)^(i) = , 

ni=l 

^((fd))-!^)^ 1 ) + fc( m )p)d«f^ m2 ) = 0. 



(175) 
(176) 

(177) 

(178) 

(179) 
(180) 



m 2 =l 



Here we combine equations f l 1 3 9 [) and (I144p into the single equation (11771) . Let us transform 
some of eqs. (ll75til80p . For instance, considering only particular solutions to eq. (j!78p we assume 



/3 (1) £o = P£o = 0. 



(181) 



In other words, the rows of (3^ and P are orthogonal to the columns of £o- Similarly, eq. (j!75j) 
means that the rows of £o are orthogonal to the columns of aS 1 ' . 

In turn, constraint (I180p may be transformed into the following one using eq. (U76|) . 



d 2 



« (1) (^ (m) - TC {m) Pa {1) )fR {m2) = 



;i82) 



m2 = l 



Since r' 1 ' is invertable, one requires M > N, which is evident due to eq. (H76j) . Then we may 
rewrite eq. (11821) without factor ct^ as follows: 



[T (m) - fC (m) Pa (1) ) fP (m2) = 0. 



(183) 



m 2 =l 



All in all, the following constraints must be satisfied: (11751 11771 



Simple examples of explicite solutions. We obtain a particular solution to the four- 
dimensional nonlinear PDE (I156p . i.e. D\ = D 2 = 2. Remember that we use double indices 

so that, for instance, (3^ = f3^ m2 \ Let M = 5, N = 2, n = 2, = = J ; £( 2 ) = 
diag(Ma), R {2) = diag(ri,r 2 ), C^ n ^ = diag(4" in2) , 4" in2) ), c[ n) = c { 2 n) = 1. To satisfy 
eqs. (11751) and (11811) we use the following structures of matrices: 



(ii) 







f(n) 



Z 

z 



Jo 

z 



z 
z 



I Z Ji 
Z J x I 



£0 



Jo 
Z 



P 



z z z z 
-j 2 z z z 









-1 











-2 1 


-1 








" K 2 


z 








Z 


Jo 




fo 




z 


z 








z 


z 








z 


z 



1 

2 



-1 1 
1 2 



;i84) 



a 



(11) 



Z 4 
K 6A 
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where Z and I are 2x2 zero and identity matrices respectively, Z4 is 4 x 4 zero matrix, 
is 6 x 4 constant matrix which will be defined below, 



J = diag(l, -1), Ji 



1 

1 



1 
-1 



Ko 



53 1 

54 S5 



(185) 



We take = diag(l, 2, 3, 4). 

To satisfy eqs. (I176j) and (j!77p (m 1 = m 2 = 1) we take 



6,1 



-(2si-3) s 2 -1 - 

1 4 

2 3 g 
-(si + 1) s 2 -1 — 

5 4 3 



si 



s 2 -2 



-si 



«2 
-«2 



--(si + 1) -s 2 

5 



1 -- 



(186) 



3 J 



r^ 11 ) = diag 0(3si-2, 3s 2 , -3. 
To satisfy eqs. (11791) and ( I183j) . m\ + m 2 > 2, we take matrices f ( niri2 ) in the following 



f (12) = 

J3 = dia; 



J 3 



f(21) 



z 



f(22) 



J 5 

z 



form: 
(187) 



-1(1 + cf% + cfVi), -(1 + cf } / 2 + cf } / 2 r 2 ) 



ri 



>'2 



J A = diag(cf 1} , -4 21) ), J 5 = diag(cf 2) , -c 2 22) ). 
Elements of matrix must be defined as follows: 

cS 12) = --(l + /,(cf 1) +cf 2 V J )), , = 1,2, 



(1* 



26 



and f^ ni ' U2 \ T^ ni ' n2 ^ are following: 

f m) = diag(i-^l,-^,ff»,ff'), 

= diag Qcf )(6 Sl + 1) + I, cf)(2 + 3 S2 ) + ff >, rf>) , 
, { 22)_ A . ( 5 + cj 21) (l + 6si)Zi 2 + cf ) (2 + 3g 2 )/ 2 . (22) „ (22 / 

r (12) = diag (-^-(5 + (1 + Qs 1 Mcf 1) + 

-1(2 + (2 + 3 S2 )/ 2 (cf ) + cf V 2 )), -r 3 (12) , -rH , 
r 2 / 

f ^ = diag f-^-(5 + (1 + Q Sl )h(cf l) + cfVx)), 
-^(2 + (2 + 3 S2 )( 2 (cf ) + cf V 2 )), -ff >, -rf > 



(189) 



kr 2 

Introduce positive parameters Pi, i — 1,2, 3, 4, by the following formulas: 

S2 =18J^T2' ^ = -5S (3Sl - 2)P " <190) 

S 4 = 3— (P3P4 - P1P2) , S 5 



(191) 



5p 3 ' 6si-4' 

so that solution V(t) reads 

hie m+m-m f 22 ( Pl e^ + p 3 e^) 
D = p x e^ + p^ 2 + ^e^ 2 + p 4 , 

/n = --(l + 6 Sl ), /12 = 2(3si _ 2) , 

/21 = ^-{vm -P 3 p 4 )(i2 Sl +p 3 - 8), / 22 = 

Here 

2 

m = vi- m, m = V2- V3, m = V4- m, v*= ^m^m^, 2 = 1,2, 3 ? 4 ? ( 192 ) 

mi,m2=l 

^-l^S.. ^.^( 4+ _JS_). (193) 



i 21 



^-(lO + / 2 cS 21) (l + 6 Sl )), 



5/ 2 

(1) _ 1 / 5 + / lC j 21) _ 104 21) (22) _ (22) S^+r.cf^ 6(ci 21) +r l C S 22) ) g ; 
° 22 ~ 5 ^ hr, +Cl iUC2 + " r 2 (4-6 Sl ) + " "" n 
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a (2) 



i {2) 

'-12 



2 + 



Ps 
6si -4' 

(l + / lC S 21) )(l + 6 Sl ) 



(194) 



+ / lC S 22) (i + 6 Sl ; 

5 + /icS 21) (l + 6 Sl 



5/ 2 (4 21) +r 2 4 zz; )(-8 +Ps + 12si) 



,(22) 



(2) _ ^ 
a 21 — ; + C 2 
f 2 



(21) 



Pa 



6s 1 - 4 



2r 2 (3 Sl -2) 



(2) 
a 22 



+ 4 21 » 



r 2 4r 2 



6r 2 si 



(3) 11 6 Sl p 3 ( 3 ) / 2 (4 21) +r 2 cf ) )(-8 + p 3 + 12gi) , 
%1 ~ T ~5~ 6^' %2 " 2r 2 (3 Sl -2) ' (195) 

4? = ^ + ^ 21)(1 + 6ai) +4 21) (2 + P3 



(3) 



a 22 



Z 2 5 ^ V 6si- 

-5/irx4 21) (-8 + p 3 + 12gi) + 2r 2 (3si - 2) (5 + /i(cj 21) + r lC j 22) )(l + 6gi)) 

10i 1 r 1 r 2 (3s 1 - 2) 



'11 



'"12 



a (4) 

a 21 



6(1 + si) 



— (5 + / 1 (cf 1) +r 1 cS 22) )(l + 6 Sl )) + 
/1 5 



(196) 



p 3 ~ / 2 (cf } + r 2 c 2 zz Q(-8 + p 3 + 12si 
2r 2 (3 Sl -2) 



(22) v 



a (4) 

u 22 



^2 



Note that not all constant matrices may be arbitrary diagonal matrices in eq. (1156j) . In fact, 
eqs. ( 1188(1 mean the following relation 



C (12) jR (2) +L (2) (C( (21) +C r( 2 2) i? ( 2 ) ) 







(197) 



Since all Pi are positive, solution V(t) (11911) has no singularities unless 



2 

mi ,m2- 



.l\t 



m\m2 I 



— > 00. However, offdiagonal elements of V tend to infinity in some directions 
in the space of parameters t mim2 . Thus, V is not bounded solution. Now we derive a simple 
example of the bounded soliton-kink solution. For this purpose we take 



S3 = s 5 = 0, s 4 = -1 
instead of eqs. (11901) . Then one gets the following formula for V: 

fll fl2 



(198) 



V(t) 



/21 /22 



(199) 



6 



11 



-(3si - 2)(6 Sl + l)s 2 , /i 2 = -3(6 Sl + l)s 2 , 
5 



hx = -2(3 Sl - 2)(3s 2 + 2), / 22 = 6(3 Sl - 2)(3s 2 + 2)s 2 , 
d = — 6(3si — 2)s 2 , 
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where rj nin2 are the linear functions of t niri2 : 
2 2 

Vl2 ^ ^ a nin2 t nirl2 , V21 ^ b nin2 t nin2J (200) 



ni,ri2=l ni,n2=l 

an = l(s x + 1), aia = — (ri(2 + (3s 2 + 2)/ 2 (cf } + cf V 2 )) + r 2 (a n - 2)), 

5 rir 2 

(21) 12 1 
fl2i = "V(6si + 1) + - - a 22 = (n(2 + (3s 2 + 2)/ 2 (cf } + cfr,)) - r 2 (a 21 / 2 + 2)), 

5 il / 2 <2^1^2 

fon = 3s 2 + 1, 6 12 = -_L(r 2 (5 + (6si + l)/i(cS 21) + c? 2) n)) + 5n(6 u - 1)), 

5rir 2 

6 21 = 4 21) (3s 2 + 2) - r + F' & 22 = T7^(r 2 (5 + (6 Sl + l)Zi(c< 21) + cf^)) - hn{b 21 h + !))■ 

h h ol 1 r 1 r 2 

If d > 0, then the diagonal elements of matrix (11991) are kinks, while offdiagonal elements tend 
to infinity in some directions in the space of parameters t mim2 . In order to obtain bounded 
solution we require 

V12 = a V2i a ij = abij, a > 0. (201) 

For the sake of simplicity, we solve eqs. (12011) for the particular choice of the arbitrary parame- 
ters: 



One has 



n =2, r 2 = 3, h =4, l 2 = 5, cS 21) = 6. (202) 

741 + 366cf ) -12cf ) (7 + 244 22) )) 
8(-559 + 108c? 2) + 36cf 2) (2 + c 2 22) ))' 
s 2 = (14521 - 1810904 22) +43200(4 22) ) 2 + 144(cf 2) ) 2 (-l + 20c 2 22) )+ (204) 

40cf 2) (-20 + 9c^ 2) + 360(4 22) ) 2 )) x 

(22)0 
2 i 



(40(-l + 15c 2 22) )(-559 + 108cj, 22) + 36cf 2) (2 + c 



246 - 36904 22) 



-233 + 2160c 2 22) + 36c5 22) (-l + 20c { 2 22) ) 



(205) 



-89 + 12cf 2) 

In additin, one has to provide positivity of a = P i > in eq. (11991) and positivity of d = p 2 > 0. 
This requirement yields constraint for q , z = 1,2: 

C J 22) = (l240 Pl - P ?(725 + 144 P2 ) - P 2 (525 + 248 P2 ) ± (207) 
Pi(l + Pi) V / 25(-2 + p 1 ) 2 -20p 1 p 2 ) (48p 2 (-5 + Pl (5 + Pa ))) ~\ 

4 22) = (- 1200 P i + 5 P 2 (-263 + 48 P2 ) + 5 P ?(497+ 100 P2 )± 
3pi(l + Pi) v / 25(-2 + Pl ) 2 - 20p lP 2 



60(10pi - 5 P ? - P 2 (15 + 2 P2 ) ± Pl (l + Pl ) v/25(-2 + Pl ) 2 -20 PlP2 ) 
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In particular, if p\ = 1/2 and pi = 1, one has 



si 



.(22) 



«2 



1 ± VT85 

24 ' 
-2545 =f 3\/l85 
~ 192" 



15 ± a/185 
30~ 



.(21) 



3(-145± 123\/185) 



5 u 2 



2080 



(208) 



.(22) 



5 "-2 



70 =f 93V185 
780 



Now expression for V, eq. (ll99|) reads 
V(f) = - 



5± 



(13 ± V185) 



20(1 + 5e( a+1 ^i) 2(e~ a ^ + 5e^\ 
(-17 ±Vl8E) (35 iv 7 !^) 



2(e-^i + 5e ar »i) 10(1 + 5e( a+1 )^i) 
Appropriate expressions for b niTl2 are following 



, f/21 = 5^ fo «m2*nin 2 - ( 2 09) 



25 , 295 ±61 v / 185 

Oil = — > °12 = — 



10 
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(210) 



3(9±2 v / l85) -214T43V/T85 

021 — — , &22 — 



10 ' 60 

We see that diagonal elements of V represent kinks, while off-diagonal elements represent non- 
symmetrically shaped solitons. 

Note that solution (I209j) causes some restrictions on the coefficients of eq. (11561) . which are 
eqs. (12061) and (12071) 



4 Conclusions 

We represent a new algorithm allowing one to construct a rich variety of particular solutions 
to a new class of nonlinear PDEs of any order in any dimensions. These equations can be 
considered as multidimensional generalizations of well known C- and S'-integrable equations. 
We show that the solution space may be rich enough to provide complete integrability of some 
of these equations. However, the problem of complete integrability requires futher study. 

The suggested algorithm allows evident generalizations. For instance, let us generalize 
constraint (150"]) as follows: 

D D 

A (m) * A {n) B {mnp) = i (m) P (mp) , p = 1, . . . , D, j = 1, . . . , mo, (211) 

m,n=l m=l 

where 5 (mrip) and P (mn) are some constant no x no matrices. Then one can show that no x uq 
matrix functions V^^t) are solutions to the following system of nonlinear PDEs: 

D D 

^2 [(v£? + (V {m) V {n) + V {m) A {n) ) ) £ (mnp) ] + ^ V [m) P {mp) = 0, (212) 

m,n=l m=l 

p=l,...,D. 

Remark that the physical application of some of the derived PDEs is obvious. For instance, 
the multidimensional iV-wave equation (J5J) appears in multiple-scale analisis of any physical 
dispersion system. However, physical applications must be considered in more details. 

The author thanks Professor P.M.Santini for usefull discussions. This work is supported by 
the RFBR grants 10-01-00787 and 09-01-92439 and by the grant NS-4887.2008.2. 
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